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Abstract 

The low infrared background and high atmospheric transparency are the prin- 
cipal advantages of Antarctic Plateau sites for astronomy. However, the poor 
seeing (between one and three arcseconds) negates much of the sensitivity im- 
provements that the Antarctic atmosphere offers, compared to mid-latitude sites 
such as Mauna Kea or Cerro Paranal. The seeing at mid-latitude sites, though 
smaller in amplitude, is dominated by turbulence at altitudes of 10 to 20 km. 
Over the Antarctic plateau, virtually no high altitude turbulence is present in 
the winter. The mean square error for an astrometric measurement with a dual- 
beam, differential astrometric interferometer in the very narrow angle regime is 
proportional to the integral of h 2 Cfj(h). Therefore, sites at which the turbu- 
lence occurs only at low altitudes offer large gains in astrometric precision. We 
show that a modest interferometer at the South Pole can achieve 10 //as differ- 
ential astrometry 300 times faster than a comparable interferometer at a good 
mid-latitude site, in median conditions. Science programs that would benefit 
from such an instrument include planet detection and orbit determination and 
astrometric observation of stars microlensed by dark matter candidates. 

Keywords: Instrumentation: Interferometers — Techniques: Interferometeric 
Astrometry — Extrasolar Planets — Galaxy: Halo, Stellar Content 



1 Introduction 

Properties of the atmosphere are the ultimate limit to astronomical observations from 
ground-based sites. Atmospheric transparency, background and seeing (indicated by 
the full width of the point spread function at half maximum, FWHM) are the param- 
eters by which sites are typically compared. For observations in the optical and near 
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infrared, high mountain sites that are above most of the Earth's boundary layer such 
as Mauna Kea and the Chilean Andes are superior because of favorable values of these 
parameters. Transparency and background are parameters that are relatively simple 
to quantify and interpret. However, there is much more to seeing than the FWHM of 
images taken through the turbulent atmosphere. 

The optical significance of turbulence in the atmosphere is due to refractive in- 
dex fluctuations, driven by temperature fluctuations. These temperature fluctuations 
are usually described in terms of a Kolmogorov model, with a temperature structure 
function as follows: 



Here, T is the temperature, x and r are location variables, and C\ characterizes the 
amplitude of the turbulence. These temperature fluctuations cause refractive index 
(N) fluctuations through the refractivity dependence of air on temperature, C% = 



The optical effect is characterized by the three-dimensional refractivity power spec- 
trum, which for Kolmogorov turbulence is 



Here, k is the spatial frequency of the turbulence, and Cfj(h) describes the vertical 
profile, where h is the height above the observatory. C%(h) is typically complicated, 
with multiple layers resulting from wind shear in the atmosphere. 

The Fried parameter, r , relates the seeing to the equations above, and is the scale 
over which the height-integrated phase fluctuation at a given wavenumber, k = 2n/\, 
equals one radian RMS. 

Given a C%{h) profile, the Fried parameter is 



and the FWHM of an image is 0.98A/r . z is the zenith angle of observation. For the 
rest of this paper, we will refer all calculations to the zenith. 

Thus, the integrated amplitude of C% determines r and therefore the seeing. The 
integrated seeing to ice level has proven to be a major drawback for Antarctic tele- 
scopes, where the seeing is one to three arcseconds at 0.5 /im, a value more typical of 
poor, low-altitude sites than superb sites such as Mauna Kea or Cerro Paranal, where 
the seeing is routinely below an arcsecond. 

However, other parameters, such as scintillation, isoplanatic angle, and tilt aniso- 
planatism, which have paramount significance in certain kinds of astronomical obser- 
vations, depend on higher order moments of the C% profile. Specifically, tilt anisopla- 
natism is the limiting error term for astrometric interferometry, and is a result of the 



D T (r) = ((T(x) - T(x + r)) 2 ) = C 2 r 2 '\ 



(SN/ST) 2 C 2 . 



$(k, h) = 0.033C 2 N (h)K 
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second moment of the C% profile. The mean square error for an astrometric measure- 
ment with a dual-beam, differential astrometric interferometer in the very narrow-angle 
regime {Oh « B) with long integrations, t » B/V, is ( |5hao fc Colavita, 1992|) : 



a] = 5.25£- 4/3 fl 2 | tfCjfWVty^h. 



In this equation, B is the baseline of the interferometer, 6 is the angular separation of 
the celestial objects, V is the wind speed as a function of height, and t is total integra- 
tion time of the observation. This calculation assumes a Kolmogorov power spectrum 
with no outer scale. Due to the h 2 factor, this integral is completely dominated by 
high altitude turbulence. Thus, a site where the Cjj profile is devoid of high altitude 
turbulence provides a substantial advantage. 

The South Pole is such a site. The C^(h) profile of the atmosphere above the South 
Pole is fundamentally different from that of any other characterised site ( [Marks et al.J 



1999| ). The unusual polar atmosphere is entirely dominated by low altitude turbulence 
and is not affected by the turbulence generating jet streams or trade winds nor the 
high-altitude synoptic wind shear that creates turbulence at mid-latitudes. 

First, we explore the benefits of a South Pole interferometer quantitatively, and 
then highlight two scientific projects for such an interferometer. 



2 Benefits of a South Pole Interferometer 

The atomspheric turbulence at a site like Mauna Kea is dominated by high altitude 
turbulence generated by tropopause instability, jet streams, and wind shear between 
synoptic weather systems. The lack of solar stratospheric heating at the South Pole 
during the polar winter results in a nearly adiabatic upper atmosphere, which means 
that the optical effects of turbulence (due to temperature fluctuations) are substantially 
reduced. The tropopause is at least a factor of two lower and very weak. Furthermore, 
there are no jet streams or trade winds at the South Pole, which are responsible for 
the large wind shears observed in atmospheric profiles at mid-latitude sites. 

In order to compare the sites, we describe the C]^(h) profile in terms of a Hufnagel- 
Valley (HV) model ([Bardy, 19981 ; IRoggemann fc Welsh, 19961) . The HV model is 



an empirically fitted heuristic model, used extensively in the context of theoretical 
calculations of light propagation through the turbulent atmosphere. The Cj^ profile 
is fitted by a series of terms: a planetary boundary layer term of the form Cfj(h) = 
A\ exp(— h/Hi), a tropopause term of the form Cjj(h) = A 2 h 10 exp(—h/H 2 ), and an 
individual layer term of the form Cjf(h) = A 3 exp( ~( h ~j? 3 ^ ). The amplitudes, A n , 
heights, H n , and the thickness, d, of the third component are all determined empirically. 
The values of these parameters used in our models of Mauna Kea and the South Pole 
are indicated in Tables |I| and ||. 
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Since the astrometric error is affected by the wind speed, we model the atmospheric 
wind profile following Greenwood (19771) , but we neglect the zenith angle and azimuth 



terms. This model consists of a constant ground layer term and a Gaussian profile 
tropopause wind term: 



V(h) = Vg round + V T exp 



h-H\ 2 



T 



For the Mauna Kea model, we use: V gronn A = 5 m s -1 , Vt = 25 m s _1 , H = 12 km, 
T = 5 km. For the South Pole model, we use a higher ground wind speed, and a more 
moderate tropopause speed at a lower altitude: V^-ound = 10 m s" 1 , Vr = 10 m s _1 , 
H = 7 km, T = 3 km. The astrometric error is not strongly sensitive to the details of 
the wind speed. 

Marks et al. (1999, |1996|) have measured the Cjy profile at the South Pole with 



tower and balloon based microthermal sensors. This profile clearly shows that the vast 
majority of the turbulence is at altitudes below 1 km. Our HV model is fit to the 
median Cff profile data of |Marks et al. (1999| ). Integrating the HV model described in 



Table |2| and Figure |2] results in seeing of 1.8 arcseconds at 0.5 fim, typical for the site 
as measured by direct imaging and differential image motion monitors. 

We now use the model to assess the astrometric errors. In order to compare sites 
and our calculations with others in the literature, we refer all calculations to an inter- 
ferometer with a 100 m baseline, reference stars separated by 1 arcminute and a 1 hour 
integration, neglecting non- atmospheric noise sources. 

For the South Pole model, the astrometric error is as = 6 /ias. We have also 
integrated a standard 0.5" seeing Mauna Kea model (Pardy, 1998|) in Table [I] and 



Figure [I]. The astrometric error for the same interferometer parameters is as = 100 /zas. 
Another commonly used reference atmosphere for Mauna Kea is that of |Koddier et al. 



(1990| ) , which is the reference atmosphere used by |Shao fc Golavita (1992| ) . Integrating 



this 0.3" seeing profile we obtain an astrometric error of as = 64 /ias, comparable to the 
Shao fc Colavita (1992j ) value of as = 60 /xas, for a 100 m baseline and a star separation 



of 1 arcminute. This mode of interferometry has been demonstrated to achieve these 
levels of precision with the Mark III and Palomar Testbed Interferometer (|GolavitaJ 
T994] ; EJolavita et al., 1999|) . 



A measure of performance is the integration time required to achieve a certain level 
of accuracy given the same interferometer configuration. Equation [I] indicates that the 
integration time is inversely proportional to erf. Therefore, these calculations indicate 
that an astrometric interferometer at the South Pole will achieve a given astrometric 
accuracy 300 times faster. A 300-fold increase in the efficiency of observation enables 
a number of science programs generally thought to be restricted to space-based inter- 
ferometers. 

This calculation is possibly optimistic. A few balloon measurements of Cfj (pVIarke 
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Table 1: Mauna Kea Atmosphere Model (|Hardy, 1998 



D- 



Term 


Boundary Layer 


Tropopause 


Strong 


Layer 


Total 




Aexp(-h/H) 


Ah 10 exp{-h/H) 


Aexp( 


-M) 2 \ 

2d? > 




A 


1 x 1(T 17 m" 2 / 3 


1.6 x 10" 53 m- 2 / 3 m- 10 


1 x 10" 


-16 m -2/3 




H (m) 


3000 


1000 


6500 






d (m) 






1000 






r (cm) 


65 


44 


27 




20 


FWHM 


0.15" 


0.23" 


0.36" 




0.5" 


as (/ias) 


23 


69 


69 




100 



The contribution of individual terms are calculated, in addition to the total. The Fried 
parameter, r is calculated at 0.5 /zm. The total r is calculated as r = (J2i r oi 5 ^ 3 ) 3 ^ 5 - 
FWHM = 0.98A/r . The calculated atmospheric astrometric error, as, is the root 
mean square error for a 100 m baseline interferometer in 1 hour on stars separated by 
1 arcminute, neglecting all non-atmospheric noise sources, as adds in quadrature. 



Table 2: South Pole Atmosphere Model in Winter Adapted from Marks et al. (1999| ) 



The terms and notation are the same as in Table |l| 
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Lower Boundary Layer 


Upper Boundary Layer 


Strong 


Layer Total 




Aexp(-h/H) 


Aexp(-h/H) 


Aexp( 


-{h-H) 2 s 
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A 


2 x 10" 15 in" 2 / 3 


2 x 10" 15 m" 2 / 3 


2 x 10" 


-14 m -2/3 


H (m) 


100 


400 


30 




d (m) 






40 




r (cm) 


21 


36 


6 


5.5 


FWHM 


0.48" 


0.27" 


1.6" 


1.8" 


as (/xas) 


2 


5 


2 


6 
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Figure 1: Mauna Kea Atmosphere Model flHardy, 1998| ). The solid line is the model 
Clj{h) turbulent atmospheric profile. The dotted line shows the cumulative contribu- 

h'^V(h')-^Cl{h')dh' 



tion to the mean square atmospheric astrometric error: 



f™h'iV(h>)~iC 2 N (h>)dh>- 
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Figure 2: South Pole Atmosphere Model in Winter Adapted from [Marks et al. (1999|) 
The terms and notation are the same as in Figure |l|. 
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et al., 1999) indicate that low- level high altitude turbulence appears at times during 
the Antarctic Winter. This turbulence is enough to degrade the relative speed to a 
factor of between 4 and 10 compared to Mauna Kea. However, this is a comparison of 
poor conditions at the South Pole to median or better at Mauna Kea. Furthermore, 
the high altitude turbulence measured by the balloon flights may be contaminated by 
wake turbulence from the balloon itself ( Marks et al., 1999| ). Although this hardly 
affects the integrated seeing, the h 2 weighting for anisoplanatism strongly emphasizes 
the high altitude turbulence. Further studies are warranted, particularly with SCIDAR 
measurements, which are most sensitive to high altitude turbulence and are immune 
to wake turbulence. 

There are two independent, though circumstantial, pieces of evidence that a large 
gain in astrometric error can be realised. First, Marks et al. ( 1999| ) calculates that 
the mean scintillation index, a], at the South Pole is four tenths of the value at Cerro 
Paranal. Since the scintillation index depends on a h 5 ^ 6 moment of the Cfj profile, we 
can estimate a lower limit to the speed improvement compared to Paranal as 0.4 -12 / 5 = 
9. 

Second, the tilt anisoplanatism term is similar to the Strehl anisoplanatism term 
that determines coherence over angular separations. Many people know the adage, 
"Stars twinkle. Planets don't." This is generally true because planets are larger than 
the isoplanatic angle at most locations on the Earth, so the scintillations of different 
parts of the planet disk are incoherent, and average out. If the isoplanatic angle 
is larger than the disk of a planet, then the planet will twinkle. During the South 
Pole winters of 1995 and 1996, one of us (JPL) saw Jupiter scintillate on occasion, 
indicating a isoplanatic angle larger than 30 arcseconds in the visible, even at an 
elevation of 22 degrees above the horizon. The mean conditions profile in Figure ^| gives 
an isoplanatic angle of 33 arcsec at 0.5 /mi at the zenith. Since the isoplanatic angle 
is proportional to airmass to the 8/5 power, this naked eye observation indicates that 
the zenith isoplanatic angle must be as large as 2 arcmin under some conditions. This 
is compelling evidence that the atmosphere above the Antarctic Plateau is uniquely 
lacking high altitude turbulence. 

The mechanism that generates the existing turbulence at the South Pole is wind 
shear across the strongly inverted lower atmosphere. Such conditions are likely to exist 
anywhere in the Antartic plateau, but the strength and height of the inversion, and 
windspeed vary markedly. The height of the inversion layer at the South Pole strongly 
correlates with wind speed, since the settling of cold air into a strong, shallow inversion 
is disrupted by high wind. Wind over the Antarctic plateau is dominated by katabatic 
flow from the high points of the plateaus, downhill towards the coast. Therefore, higher 
sites such as Dome C, where the mean windspeed is only 2.6 m s _1 are likely to show 
strong, shallow thermal inversions, confining the turbulence generated to even lower 
altitudes, which will offer even larger gains to astrometric interferometry. Likely to be 
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the ulitimate site is the highest point on the plateau, Dome A. At these sites it is even 
possible that the turbulent inverted layer may be shallow enough that it is feasible to 
elevate the telescope above it, and realise the "super seeing" suggested by |Gillinghamj 



(1993Q . Site testing programs are underway to characterise these sites, including seeing 
monitoring programs. Measurements of the Cff profiles will be crucial to determination 
of the scientific programs and instruments to be deployed at these sites. 



3 Science Applications 

The benefits shown in this calculation are limited to differential astrometric obser- 
vations. Such observations have many applications. We highlight two of them here, 
which are particularly promising with regard to important astrophysical questions and 
well suited to be conducted in the southern hemisphere. We have also considered other 
configurations to take advantage of the unique properties of the Cjf(h) profile. 

The factor of 300 in speed achievable with an interferometer at the South Pole 
result from the low altitude of the turbulence causing common aberrations that can 
be removed with a differential measurement. The same technique can be exploited 
in different ways with single telescopes and adaptive optics or interferometers. The 
large isoplanatic angle at the South Pole is attractive for adaptive optics applications. 
However, the small ro requires a very bright guide stars and a high order AO system, 
even for a moderate sized telescope. Similarly, a conventional interferometer measuring 
fringe visibility suffers from the poor seeing. A dual star interferometer is required to 
gain an advantage from the low altitude of the turbulence with phase referencing, 
taking advantage of the coherence of the atmosphere across large angles. Such an 
interferometer can be used to measure fringe visibility on a faint source, while tracking 
on a bright source. A generally useful interferometer in this mode would require a large 
number of telescopes and long baselines to allow bootstrapping the phase referencing 
to long baselines. To measure relative fringe phase, however, only requires a single 
baseline, and it is disadvantageous to resolve the source, so moderate baselines are 
desirable. We therefore conclude that astrometric interferometry is uniquely suited to 
the site characteristics and logistical constraints in Antarctica. 



3.1 Planet Detection and Orbit Determination 

Astrometric observations with accuracies near a few //as permit the detection and 
determination of the orbits of extrasolar planets. For reference, if a star like the Sun, 
lying at a distance of 10 pc, has a Jupiter-like planet orbiting it, the star's position 
will be perturbed by a maximum of 1 mas over the orbital period (~ 12 years). The 
signature of an Earth-like planet is approximately 1 //as. Because the South Pole 
interferometer could, in principal, conduct observations indefinitely, within 10 years 



9 



one could not only survey and determine the three dimensional orbital characteristics 
of the planets discovered through radial velocity studies, but also survey a substantial 
number of other stars within 30 pc for signatures of extrasolar planets. By comparison, 
the Space Interferometry Mission (SIM), which has a limited lifetime of 5 years, will 
be unable to observe stars over such long periods, limiting the types of planets it will 
find. It is necessary to fit the mass, semi-major axis, eccentricity, parallax and proper 
motion of the star in the astrometric solution. Accurate determination of all of these 
parameters with less than a complete orbit is not possible. It will therefore require 
significantly longer temporal baselines than SIM can provide to fully characterise the 
distribution of planets. 

In comparison to the Keck Interferometer, the advantages of the South Pole inter- 
ferometer are the factor of 300 in speed of observation, meaning that many more stars 
could be surveyed in a given time. Alternatively, the substantially larger isoplanatic 
angle permits observation of more stars. For example, instead of speed, the benefits 
result in a 300-fold increase in the area of sky from which reference stars can be chosen. 



3.2 Astrometric Measurement of OGLE Microlensing Events 

There is a controversial body of evidence that suggests that a substantial fraction of 
the dark matter in the Milky Way is in the form of old white dwarfs. These may have 
been detected by microlensing experiments (MACHO; |Alcock et al. (2000| )) and direct 
surveys ( pppenheimer et al., 200"T| ). However, the locations of the lensing masses re- 



mains unmeasurable and is a major source of contention among various researchers. It 
is possible that none of the microlensing events correspond to lenses in the Galaxy's 
halo. Astrometry of the source star during a microlensing event breaks the degeneracy 
in the solution to the microlensing equations. As a result, one can determine the mass 
and distance to the microlens itself. As shown by |Boden et al. (1998|) , astrometric 



accuracy of better than 10 /xas permits accurate determination of the lens parameters, 
although some degeneracy remains concerning the transverse velocity of the lens. In- 
deed, a lens placed at a distance of 8 kpc with the source star in the Large Magellanic 
Cloud (LMC), produces an astrometric effect of about 100 /xas, readily detectable with 
a modest baseline South Pole interferometer. 

Although the MACHO project QAlcock et al., 200q ) is no longer searching for mi- 



crolenses, the OGLE-II project ( [Udalski et al., 1997|) is currently engaged in more 
sensitive observations of the LMC and their system for broadcasting early detection of 
new events is almost operational. The LMC is easily observed from the South Pole. 
With continuous observation possible during the polar night, an interferometer at the 
South Pole could resolve the remaining question concerning the nature of the MACHO 
objects. 
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4 Conclusions 



We have demonstrated that the unique structure of the atmosphere above the high 
Antarctic plateau provides significant advantages to interferometric astrometry. The 
calculations described above motivate improved measurements of the Cfj profile of 
the atmosphere with statistical coverage of varying conditions. If such measurements 
confirm our results, the construction of an interferometer at the South Pole should be 
considered for the near future. The observations by such an instrument will have a 
substantial impact on several important astrophysical problems. 
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